Strongly driven molecules: probing the tunneling phase in the over-the-barrier regime 
and prevalence of different double ionization pathways 
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Using a three-dimensional quasiclassical technique we explore molecular double ionization by a 
linearly polarized, infrared (800 nm) 27 fs laser pulse. For intensities ranging from the tunneling to 
the over-the-barrier regime, we identify the double ionization pathways in a unified way as a function 
of total electron energy. For the tunneling regime, we discuss the differences in the interplay of double 
ionization (DI) pathways between strongly driven He and strongly driven N2. For intermediate 
intensities in the over-the-barrier regime, we find that both the correlated momenta and the double 
ionization probability distribution as a function of total energy probe the tunneling phase of the 
re-colliding electron. This allows for a direct verification of the re-colliding electron tunneling at a 
large phase of the laser field in the over-the-barrier regime in contrast to a small tunneling phase in 
the tunneling regime. 



PACS numbers: 33.80.Rv, 34.80.Gs, 42.50.Hz 



I. INTRODUCTION 

Over the past two decades electron correlation has 
been established as the tinderlying mechanism for many 
important phenomena arising from the interaction of 
strong laser pulses with matter. One of these phenomena 
is the dramatically enhanced multiple ionization yield of 
atoms (e.g. [l|) and molecules (e.g. (H) in intense laser 
pulses. When the field intensity is large the two electrons 
are stripped out sequentially (3| while for smaller intensi- 
ties non-sequential double ionization (NSDI) dominates. 

According to the accepted mechanism for NSDI — the 
three-step or re-scattering model [1] — 1) one electron es- 
capes through the field-lowered Coulomb potential, 2) it 
moves in the strong infrared laser field and 3) it returns 
to the core (possibly multiple times) to transfer energy 
to the other electron in the parent atomic or molecu- 
lar ion. Using coincidence imaging techniques such as 
COLTRIMS many experiments have succeeded in ob- 
taining highly differential kinematical details of electron 
correlation in the non-sequential intensity regime (e.g. 
[^-[U). However, for higher intensities the available ex- 
perimental data is much less abundant. Correspondingly, 
the majority of theoretical work has concentrated on the 
non-sequential intensity regime (e.g. [1243 )• 

Addressing the double ionization of strongly driven 
systems with first-principle quantum mechanical calcu- 
lations in all three spatial dimensions (3-d) is an im- 
mense task. Currently, 3-d ab-initio quantum mechan- 
ical calculations are mostly available for the driven He 
15| . With the exception of few 3-d ab-initio stud- 
io, to cope 
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with the highly complex task of tackling double ioniza- 
tion of molecules, many studies use numerical quantum 
approaches of reduced dimensionality. Others use judi- 
ciously chosen quantum models of reduced dimensional- 
ity (e.g. [l7|). while some use semi-analytical quantum 



approaches in the framework of the so-called Strong-Field 
Approximation (SFA) (e.g. (3), not fully accounting for 
the Coulomb singularity. 

Here, we report on a classical study of electron correla- 
tion in the N2 molecule driven by a 27 fs laser pulse at 800 
nm for intensities well within the non-sequential double 
ionization regime as well as for intensities corresponding 
to the transition from pure tunneling to over-the-barrier 
ionization. We use a three-dimensional quasiclassical 
technique that we developed for conservative systems, i.e. 
the single photon multiple ionization of atomic systems 
[T^-Hlj . We extended this technique to non-conservative 
systems to first treat the correlated electron dynamics of 
the strongly driven He atom 



strongly driven N2 molecule 
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and subsequently of the 
. Our method is numer- 



ically very efficient and treats the Coulomb singularity 
with no approximation in contrast to techniques that use 
"soft-core" potentials. 

We describe in a unified way the different double ion- 
ization (DI) pathways the two electrons follow to escape 
after re-collision as a function of total electron energy for 
the tunneling regime and for intensities corresponding 
to the transition from pure tunneling to over-the-barrier 
ionization — we refer to the latter as intermediate inten- 
sities. (Throughout this work the total electron energy 
we refer to is the sum of the final kinetic energies of the 
two electrons). A similar description was very recently 
used to describe the DI pathways of strongly driven He 
[U, [23 . This general treatment allows us to identify uni- 
versal features of strongly-driven molecular double ion- 
ization and compare them with the atomic case. In a 
previous study on strongly driven N2 [23| . we identified 
the DI pathways only in the tunneling regime and inde- 
pendent of total electron energy for a very short 7 fs laser 
pulse. In our current study, we drive N2 with a long laser 
pulse. 

We first discuss the general features of DI for differ- 
ent total energies and for intensities in the tunneling 
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and the over-the-barrier regime. We then show that, for 
all intensities considered, the DI pathways have distinct 
patterns of correlated momenta along the field polariza- 
tion axis. Interestingly, for intermediate intensities in the 
over-the-barrier regime the correlated momenta of the Di- 
rect pathway — simultaneous ejection of both electrons — 
have an anticorrclation pattern. Moreover, the correlated 
momenta of the Delayed pathway — ejection of one elec- 
tron with a delay of more that a quarter of a cycle with 
respect to re-collision — have a square pattern. Both the 
anticorrelation and square patterns arise from the "soft" 
collisions, re-collision assisted field ionization [2^, under- 
lying the intermediate intensities in the over-the-barrier 
ionization. We show that the electron momentum corre- 
sponding to the boundary of the square pattern probes 
the tunneling phase of the re-colliding electron. 

Finally, we explore in detail the features of the DI prob- 
ability distribution as a function of total electron energy. 
We find that in the tunneling regime the DI probability 
distribution has one peak, while for intermediate intensi- 
ties in the over-the-barrier regime it has two peaks. We 
show that one peak — present for all intensities — is due 
to the field-assisted re-collision mechanism, and shifts 
to smaller energies with increasing intensity. The other 
peak is due to the "soft" collisions underlying DI for 
intermediate intensities in the over-thc-barrier regime. 
This latter peak probes in yet another way the tunneling 
phase. 



II. METHOD 

The method we use was previously described in [23| . 
For completeness we also describe it in what follows. Our 
3-d quasiclassical model entails the following steps: We 
first set-up the initial phase space distribution of the two 
"active" electrons in the N2 diatomic molecule. Here, 
we consider only parallel alignment between the molec- 
ular axis and the laser electric field. At intensities in 
the tunneling regime we assume that one electron tun- 
nels through the field-lowered Coulomb potential. For 
the tunneling rate one can use quantum mechanical or 
semiclassical formulas for diatomic molecules (see, e.g. 
[H-iil). We use the rate provided in ref. The 
longitudinal momentum is zero while the transverse one 
is provided by a Gaussian distribution [soj . This de- 
scription is valid as long as the potential barrier is not 
completely suppressed by the instantaneous laser field 
E{t) = Eo{t) cos{Lut). We consider the usual laser wave- 
length of 800 nm, corresponding to w = 0.057a.u. (a.u. 
- atomic units). In our simulation the pulse envelope 
Eo{t) is defined as Eo{t) = Eq toi < t < 7T and 
Eo{t) = Eo cos2(w(t - 7r)/12) for 7T < t < lOT with T 
the period of the field. The threshold for over-the-barrier 
ionization in neutral A^2, with an ionization energy of 
/pi = 0.5728 a.u., is reached at a field strength oi E = 
0.075a.u. (corresponding to roughly 2 x lO^^'Watts/cm^). 

Above 2 X 10^^ Watts/cm^ the laser field allows an un- 



hindered electron escape and therefore the initial phase 
space is modeled by a double electron microcanonical dis- 
tribution [U. However, in setting-up the initial phase 
space distribution we transition smoothly from the tun- 
neling to the over-the-barrier intensity regime. Namely, 
we assign a random number to the phase (j) of the laser 
field when the first electron is ionized, see [H, [s^. If the 
phase corresponds to an instantaneous strength of the 
laser field E{(j)) that leaves the electron below the barrier 
then we use the initial conditions dictated by the tunnel- 
ing model. If the instantaneous field strength pushes the 
barrier below the Ipi of that electron then we use the mi- 
crocanonical distribution to set-up the initial phase space 
distribution. This choice of initial conditions has proven 
successful in past studies |33j in modeling the experimen- 
tal ratio of double versus single ionization for long laser 
pulses [34} . With our approach we ensure a smooth tran- 
sition of the initial phase space distribution as we change 
the intensity. Even at an intensity of 3x lO^^Watts/cm^ 
still about 70% of the double ionization probability cor- 
responds to trajectories initialized using the tunneling 
model, while 30% of the probability corresponds to tra- 
jectories initialized using the microcanonical distribution. 

After setting-up the initial phase space distribution we 
transform to a new system of coordinates, the so called 
"regularized" coordinates This transformation is 

exact and explicitly eliminates the Coulomb singular- 
ity. This step is more challenging for molecular systems 
since one has to "regularize" with respect to more than 
one atomic centers versus one atomic center for atoms. 
We regularize using the global regularization scheme de- 
scribed in ref. [s^. Finally, we use the Classical Trajec- 
tory Monte Carlo (CTMC) method for the time propa- 
gation [13. The propagation involves the 3-d four-body 
Hamiltonian in the laser field with "frozen" nuclei: 
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where E{t) is the laser electric field polarized along the 
z direction and further defined as detailed above, and R 
is the internuclear distance. 



III. RESULTS AND DISCUSSION 

We consider three laser intensities at 10^^, 1.5 xlO^^ 
and 3 xlO" Watts/cm^. The range of laser intensi- 
ties we consider is important because, at 800 nm, for 
1.44x lO^'' Watts/cm^ the maximum return energy of the 
re-colliding electron (3.2 Up according to the three-step 
model equals the first ionization energy of the ground 
state of ■ The ponderomotive energy Up — Eq/{4:Uj'^) 
is the cycle-averaged energy of the oscillatory motion. 
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FIG. 1. We plot (pi.z) (thick black line), {p2,z) (thick grey 
line), (-Z/\R/2-ri\- Z/\-R/2-ri\) (thin black line) 



a.nd(^~Z/\R/2~r2\~ Z/\- R/2 - r2|^(thin grey line). The 

time axis is measured in laser cycles. Electron 1 is the re- 
colliding electron, fi is the position vector of electron i, and 
Z = 1. The plotted quantities are shown for 1.5 xlO^'* 
Watts/cm^: top row for energies [0,1] f/p, bottom row for en- 
ergies [8,9] (7p. Left column refers to the SE pathway and right 
column to the RESIb pathway. 
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FIG. 2. First and third columns from the left same as Fig. [T] 
but for 3 xlO^"' Watts/cm^: top row for energies [0,0.4]L'j,, 
bottom row for energies [3.2, 3.6]Up. We also show the proba- 
bility distribution for the phase when the re-colliding electron 
tunnels out for the SE pathway, second column from the left, 
and for the RESIb pathway, fourth column from the left. 



A. Mean values 

In this section, we briefly explore the characteristics 
and general properties of the double ionization (DI) path- 
ways. Our goal is to point out that for different intensities 
and different total energy regimes the DI pathways have 
different characteristic features. 

To identify the main DI energy transfer pathways we 



use the time delay between the re-collision time and the 
time of ionization of each electron [2^. We define the 
re-collision time as the time of minimum approach of the 
two electrons. We identify this time through the maxi- 
mum in the electron pair potential energy. The ionization 
time for each electron is defined as the time when the sum 
of the electron's kinetic energy (using the canonical mo- 
mentum) and potential energy with the two atomic cen- 
ters becomes positive and remains positive thereafter — 
for more details for the time of ionization see [22| and 
references there in. 

Following the steps outlined above, we identify the fol- 
lowing pathways: the Direct (SE) and Delayed (RESI) 
which are well established. (The delayed pathway is also 
referred to as re-collision-induced excitation with sub- 
sequent field ionization, RESI [H, H^). In the direct 
ionization pathway (SE) both electrons are ionized si- 
multaneously very close (less than a quarter of a laser 
period) to the rc-collision time. In the delayed ionization 
pathway, the rc-coUiding electron excites the remaining 
electron but does not ionize it. The electron is subse- 
quently ionized at a peak (RESIa) or at a zero (RESIb) 
of the laser electric field. In addition to the SE and RESI 
pathways we also identify the double delayed pathway 
(DDE) with both electrons ionizing more than a quarter 
of a laser cycle after re-collision. This pathway was first 
identified in our recent study of strongly driven He [25|. 

In Fig.[T]for 1.5x lO" Watts/cm^ and in Fig.[2]for 3x 
10" Watts/cm^, we plot the average of the momentum 
component of each electron along the polarization axis 
(the z-axis), < pi,z > and < p2,z >, for the SE and 
RESIb pathways. In addition, we plot the average of the 
potential energy of each electron with the two nuclei. The 
time of re-collision is identified in Fig. [1] and Fig. [2] as the 
time at which there is a sudden rise/dip of the nuclear 
contribution to the potential energy of the second/first 
electron. In the SE pathway, electron 2 ionizes at a time 
close to the time of re-collision while in RESIb electron 
2 ionizes around T/2 later. In both Fig. [T] and Fig. [2] 
we consider only those trajectories where the re-collision 
occurs at the first return of the re-colliding electron to the 
atomic centers. For 800 nm, we find this to be the most 
important contribution to the SE and RESI pathways for 
all intensities. Nevertheless, multiple returns of the re- 
colliding electron are explicitly accounted for in the rest 
of our results. In addition, the general properties of SE 
and RESIb pathways for one return of the re-colliding 
electron, described below, hold true for multiple returns 
the only difference being the re-collision time. 

Tunneling regime. For intensities in the tunneling 
regime above 1.44x 10^^ Watts/cm^, we find that the 
re-collision time shifts from T/2 for small total energy 
to (2/3) T (the time of maximum energy re-collision in 
the three-step model) for large total energy — see Fig. [T] 
and compare the re-collision times of a) and b). This is 
similar to our finding for strongly driven He for intensi- 
ties where 3.2 Up is above the first excitation energy of 
He"*" [1^. The behavior of < pi^z > against time also 
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differs markedly as we go from one energy regime to the 
next. For small total energy, at the instant of re-collision 
the momentum of the re-colliding electron first decreases 
due to transfer of energy to electron 2 — dip in < pi^z > 
around T/2, see the SE pathway depicted on the top row 
of Fig.[TJ But at this early re-collision time the laser field 
docs not undergo a sign change and so continues to pull 
electron 1 in the same direction as prior to rc-collision. 
For high total energy, the interaction of the rc-colliding 
electron with the nucleus becomes very strong. At 2/3T 
(re-collision time at high total energies) electron 1 gets 
sharply pulled back both by the laser field and the nu- 
cleus and reverses in direction, see bottom row of Fig. [T] 
Over-the-barrier regime. For intermediate intensities 
in the over-the-barrier regime, for the majority of tra- 
jectories, re-collision takes place early at T/2, see Fig. [2] 
Moreover, as shown in Fig. [21 we find that the re-colliding 
electron for the majority of DI trajectories tunnels when 
the phase of the laser field is around 50° — much larger 
than 17°, the phase that according to the three-step- 
model corresponds to maximum return energy. At 3x 
10^^ Watts/cm^, tunneling at a large phase of the field 
around 50° results in a maximum excursion of electron 1 
in the laser field, before its return to the molecular ion, 
which is smaller than the maximum excursion at 1.5 x 
10^^ Watts/cm^ — small tunneling phase in the tunneling 
regime. (That larger tunneling phase results in a smaller 
maximum excursion can be easily seen using the simplest 
version of the three-step model.) Indeed, focusing on the 
potential energy of electron 1 in Fig. [2] and Fig. [U we 
find that at T/4 the potential energy is smaller at 3x 
10" Watts/cm^ than at 1.5x lO" Watts/cm^. More- 
over, upon its return to the molecular ion, electron 1 
transfers a smaller amount of energy to electron 2 at 3 x 
10" Watts/cm^ than at 1.5 x lO" Watts/cm^. Indeed, 
the change in the potential energy of electron 1 during 
re-collision is smaller at 3x 10"'^^ Watts/cm^ than at 1.5 x 
10" Watts/cm^, compare Fig.[2]and Fig.[Tl The fact that 
the two electrons undergo "soft" collisions at 3x 10^^ 
Watts/cm^ is further supported by our finding that the 
re-colliding electron ionizes before the re-collision time 
for 30% of the SE trajectories and for 70% of the RESI 
ones. Thus, in contrast to "soft" collisions underlying DI 
for intermediate intensities in the over-the-barrier regime, 
it is strong collisions with both electrons ionizing after re- 
collision that underly DI in the tunneling regime. 



B. Correlated Momenta 

In this section, we further explore the presence of 
"soft" collisions at an intensity 3 x lO^''^ Watts/cm^ and 
identify their traces on the correlated momenta for each 
DI pathway. 

In Fig. [3] we show the correlated momenta of the two 
electrons for the SE, RESIa. RESIb pathways as well as 
the total pattern with all DI pathways combined. We 
find that the RESIa and RESIb pathways have distinct 



patterns only for 10^^ and 1.5 x 10"'^^ Watts/cm^ while for 
3 X 10^^ Watts/cm^ they are almost indistinguishable. 
For intensities in the tunneling regime the time electron 2 
ionizes is around an extremum of the field for RESIa and 
a zero of the field for RESIb. However, in the over-the- 
barrier regime electron 2, in the Delayed pathway, ionizes 
mostly half way between an extremum and a zero of the 
field. Clearly, at the latter time the intense field lowers 
sufficiently the Coulomb potential to still allow electron 
2 to escape. Thus, in the following sections, we only 
present the combined contribution of RESIa and RESIb 
at 3 X 10" Watts/cm^. 



We further note that at 3 x lO" Watts/cm^ the corre- 
lated momenta of RESI have a pronounced square struc- 
ture, see Fig. [21 This feature can be understood in terms 
of the "soft" collisions that underly DI. For 70% of the 
RESI trajectories the re-colliding electron ionizes at a 
very early time, before re-collision at T/2. For the latter 
trajectories the final momentum of the re-colliding elec- 
tron is mostly determined by the vector potential at the 
time the electron tunnels out (first step in the three-step 
model). At 3 x lO" Watts/cm^, elect ron 1 tunnels out at 
phases of the field centered around 50° and extending up 
to 60°, see Fig. [2l The momentum an electron gains by 
the field when the initial phase is around 50-60° is consis- 
tent with the value of the momentum at the boundary of 
the square pattern of the RESI pathway. Since the RESI 
pathway prevails overall for strongly driven N2, see fol- 
lowing section, it is the square pattern that prevails in 
the combined correlated momenta, see Fig. [31 Thus, the 
correlated momenta probe the tunneling phase of the re- 
colliding electron for intermediate intensities in the over- 
the-barrier regime. This is not the case in the tunneling 
regime; the two electrons undergo strong collisions and 
the final momenta for both electrons are determined by 
the vector potential at the re-collision time. 



The SE correlated momenta have a pronounced anti- 
corellation pattern as shown in Fig. [31 For the majority 
of SE trajectories the re-colliding electron tunnels at a 
phase around 50°. For 30% of the SE trajectories the 
re-colliding electron ionizes at an early time before re- 
collision at T/2. The second electron ionizes just after 
re-collision at T/2. Thus the two electrons ionize at times 
corresponding to opposite signs of the vector potential 
resulting in opposite final electron momenta. For the re- 
maining 70% of SE events both electrons ionize mostly 
after T/2. Out of these 70%, for some SE events the 
two electrons escape in the same direction, as is the case 
in the tunneling regime. However, for others, a small 
amount of energy transfer and possibly the Coulomb re- 
pulsion, result in the two electrons escaping in opposite 
directions. Overall electrons escaping in opposite direc- 
tions is the biggest contribution to SE giving rise to an 
anticorrclation pattern. 
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FIG. 3. Correlated momenta at intensities 10^* Watts/cm^ 
(top row), 1.5 X 10^* Watts/cm^ (middle row), and 3 x lO" 
Watts/cm'^ (bottom row). Columns from left to right indi- 
cate the various DI mechanisms: SE, RESIa, RESIb, and the 
contribution of all DI mechanisms combined. 



C. Double ionization probability distribution 

In this section we show that for 3 x 10^"' Watts/cm^, in 
addition to the correlated momenta, the DI probability 
distribution as a function of total electron energy probes 
the tunneling phase of the re-colliding electron. 

a) % contribution of DI pathways. In Fig. |4l top row, 
we plot the probability distribution of DI for 10"'^"', 1.5 x 
10^'^ and 3x 10^'^ Watts/cm^. Moreover, for each inten- 
sity, we plot the % contribution of each DI pathway as a 
function of the total electron energy, see bottom row in 
Fig. 01 We use energy steps as small as the immense com- 
putational challenge of the endeavor allows — 1-2 million 
DI events for the whole energy regime. 

Small total energies. For small intensities the DDE 
pathway prevails for small total electron energy with a 
% contribution changing from 80% for 10^^ Watts/cm^ 
to 30% for 3.0 X lO" Watts/cm^. We find that for N2 in 
the DDE pathway the re-colliding electron does not get 
significantly trapped by the core after re-collision and 
thus the two electrons do not escape in opposite direc- 
tions as is the case for strongly driven He, see [2^. A 
possible explanation is that the presence of two nuclear 
centers instead of one with nuclear charge of one instead 
of two (He charge) makes it harder for the re-colliding 
electron to become bound after re-collision. As the in- 
tensity increases to 3.0 x 10^** Watts/cm^ crossing over 
to the over-the-barrier regime, for small total energy, the 
sequential ionization (SI) pathway takes over. We define 
SI as the pathway where the two electrons are closest to 
each other at the beginning of the propagation. Thus the 
SI pathway is a signature of sequential ionization with 
both electrons ionizing mainly due to their interaction 
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FIG. 4. Top row: Double ionization probability distribution; 
Bottom row: % contributions of the main DI mechanisms (• 
for DDE, -4 for RESIa, ■ for RESIb, * for SE). First column: 
10" Watts/cm^ second column: 1.5 x lO" Watts/cm^, third 
column: 3 x 10^'' Watts/cm^ (at this intensity an additional 
DI mechanism appears, the SI mechanism, given by ♦). Note 
that for 3 x 10^* Watts/cm'^ we only show the combined con- 
tribution of RESIa and RESIb. 

Intermediate and high total energies. For all intensities 
considered, the Delayed pathway prevails for all energies, 
with the exception of very small ones. The % contribu- 
tion of SE is always smaller than that of the RESFs. A 
trace of collisional physics on strong field ionization is 
the overall shape of the SE contribution to DI as a func- 
tion of total energy. It is very small for small and high 
energies while it is large for intermediate energies, as we 
have already found for strongly driven He [25|. Thus, SE 
qualitatively resembles electron impact ionization. 

For strongly driven He, RESI prevails independent of 
the total electron energy [2J1 for intensities where 3.2 Up 
is below the first excitation energy of He"*". For higher 
intensities, for strongly driven He, the SE pathway pre- 
vails for intermediate total electron energy, see 12J, [23 . 
In contrast, for strongly driven iV2, RESI prevails for all 
intensities. The larger contribution of RESI to strongly 
driven 7V2, possibly, arises because the interaction of the 
field with the molecular ion makes available a much larger 
number of excited states to the bound electron compared 
to the atomic case. Upon return of the re-colliding elec- 
tron to the molecular ion, electron 2 can be found in any 
of these excited states. It is thus less probable for elec- 
tron 2 to be found close to either of the nuclei and to 
interact strongly with electron 1, escaping through SE; 
it is rather more probable for electron 2 to gain some 
energy from the re-colliding electron, get further excited, 
and subsequently ionize. 

We also note, that the % contribution of different DI 
pathways depends on the duration of the pulse. Our find- 
ing in the current work that RESI prevails for strongly 
driven N2 stems from the fact that we are driving N2 
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FIG. 5. Double differential probability in energy: top 
row for 10^* Watts/cm^ for the energy intervals [0,l]Up, 
[l,2]Up, [2,3]Up, [3,4]Up and [4,5]Up; bottom row for 3 x lO" 



Watts/cm^ for [0,0.8]Up, [0.8,1.6]Up, 
and [3.2,4]Up. 



[1.6,2.4]Up, [2.4,3.2]Up 



with a long laser pulse as is the case in our recent studies 
of strongly driven He [ll, [H] . 

b) DI probability distribution in the over-the-barrier 
regime. For 3x 10"'^"' Watts/cm^, we find that the DI 
probability distribution has two features that are not 
present for intensities in the tunneling regime. First the 
DI probability distribution has two peaks. Second the 
probability for two electrons to ionize with a large to- 
tal energy is very small. Fig. 2] shows that for 3x 10^^ 
Watts/cm^ the probability for two electrons to escape 
with a total energy of 6 Up is comparable to the proba- 
bility for two electrons to escape with a total energy of 
lOUp when the intensity is 10"^^ or 1.5x 10^"* Watts/cm^. 
In what follows we explain the above two features. 

As the intensity increases in the tunneling regime the 
peak of the DI probability distribution shifts to smaller 
total electron energy, from 2.5 Up at 10^"* Watts/cm^ 
to 1.6 Up for 1.5x10" Watts/cm^, see Fig. HI We also 
find that around the maximum of the DI probability dis- 
tribution the two electrons share the energy in all pos- 
sible ways, see top row of Fig. [S] The electrons share 
the energy in a similar fashion for 3x 10^^ Watts/cm^ 
when the total electron energy is around 0.7 Up — the lo- 
cation of the first peak in the DI probability distribution. 
Thus, the peak that is present for all laser intensities 
considered shifts to smaller total energy with increasing 
intensity. This is consistent with a less efficient energy 
transfer from electron 1 to electron 2, case in point the 
"soft" collisions that underly double ionization for the 
larger intensity considered in this study. The same rea- 
soning holds for the overall shift of the DI probability 
distribution to smaller total energy. 

We next explore the origin of the second peak around 
1.9 Up for 3x 10" Watts/cm^. From Fig. [5] (bottom 
row) we find that, when the total energy is between 1.6- 
2.4 Up, one electron — we find it to be the re-colliding 
one — escapes with most of the energy while electron 2 
escapes with a very small part of it. This asymmetric 
energy sharing implies that the second peak in the total 
DI probability distribution is related to the re-colliding 



electron's probability distribution as a function of its fi- 
nal energy. We check whether this is the case by plot- 
ting, in Fig. [21 for each DI pathway the probability for 
each electron to escape with a certain amount of energy. 
We find that it is indeed the first electron's RESI prob- 
ability distribution that accounts for the peak at higher 
total energy. Namely, in Fig. [5] we show that the RESI 
probability distribution of electron 1 , as a function of the 
electron's final energy, varies smoothly for intensities in 
the tunneling regime (bottom row), in contrast, to its 
pronounced peak around 1.6 Up for 3 x lO^'* Watts/cm^ 
(top row) . To further understand the origin of this peak 
we also explore the dependence on total energy of the 
70% of RESI events where electron 1 ionizes very early 
on, before re-collision at T/2. We find that the % contri- 
bution of the latter trajectories increases monotonically 
from 8% at 0.4 Up total electron energy to its maximum 
value of 80% for total energy at and above 1.6 Up. We 
conclude that this second peak in the total DI probabil- 
ity distribution is directly related to these 70% of RESI 
events. Indeed, the tunneling phase of 50-60° results in 
a final kinetic energy of 1.5 Up very close to the peak of 
the re-colliding electron's RESI probability distribution 
shown in Fig. [6l Thus this second peak in the DI prob- 
ability distribution is yet another probe of the tunneling 
phase. 




0.0001 



energy (Up) 



FIG. 6. Top row 3 x lO" Watts/cm^ first column from the 
left: total DI probability distribution as a function of total 
energy (circles) and probability distribution as a function of 
one electron's final energy, for the re-colliding electron (black 
line) and for the second electron (grey line) . Second and third 
column same as the first but for the SE and RESI pathways. 
Bottom row, same as top row but for 10^^ Watts/cm^. 



IV. CONCLUSIONS 

We have explored in detail the prevalence of the DI 
pathways for strongly driven N2 for intensities ranging 
from the tunneling regime to intermediate ones in the 
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ovcr-the-barrier regime. We have found that the Delayed 
pathway prevails for all intensities and for all total en- 
ergies, with the exception of very small energies. This 
is different from the atomic case where SE prevails for 
intermediate total energies for intensities where 3.2 Up is 
above the first excitation energy of Hc+ . 

For intermediate intensities in the over-the-barrier 
regime, wc find that "soft" collisions undcrly DI. In these 
"soft" collisions electron 1 tunnels at a large phase of the 
laser field, around 50°, has a small maximum excursion 



in the laser field, and upon its return to the molecular 
ion transfers a small amount of energy to electron 2. One 
feature of "soft" collisions is that electron 1 ionizes early 
on, before re-collision at T/2 — this holds true for 70% of 
RESI and for 30% of SE trajectories. For the interme- 
diate intensities in the over-the-barrier regime, we have 
shown that the boundaries of the square pattern of the 
correlated momenta as well as one of the two peaks of the 
DI probability distribution as a function of total energy 
probe the tunneling phase of the re-colliding electron. 
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